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Abstract 



We estimate the survival probability of J /ip in relativistic heavy-ion col- 
lisions, using both short-distance QCD and nuclear absorption mechanism. 
The suppression is found to be almost 100 percent for a thermally equili- 
brating quark-gluon plasma. The measurement of such a huge suppression of 
J /if) may suggest the existence of a deconfined partonic medium, possibly a 
thermalised quark-gluon plasma, at RHIC. 
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I. INTRODUCTION 



Recently a lot of effort is being made to detect a new phase of matter, namely quark- 
gluon-plasma(QGP). It is known from lattice quantum chromodynamics |I| that hadronic 
matter at sufficiently high temperatures (~ 200MeV) and densities undergoes a transition 
to this quark-gluon plasma phase. While such a phase surely did exist in the early universe, 
it is interesting if we can recreate the early universe and experimentally verify this QCD 
phase in the laboratory i.e. in ultra relativistic heavy ion collisions (URHIC). In the near 
future, the relativistic heavy ion colliders(RHIC) at BNL and large Hadron Colliders(LHC) 
at CERN will be able to study such a phase. However this deconfined phase is not 
accesible directly at these collider experiments. The signatures are therefore necessarily 
indirect, and the prominent ones are : 1) J/ip suppression Q, 2) electromagnetic probes 
such as dilepton and direct photon production and 3) strangeness enhancement Q]. 

Dileptons and single photons have long been proposed as useful probes of the plasma [0 , 
as once produced, they hardly interact with the strongly interacting matter and thus carry 
the details of the circumstances of their production. However, as they are also produced 
via hadronic decays in the later stages of evolution, where plasma expands and cools, it 
becomes difficult to distinguish between them. In that sense the production of cc is a clean 
process, as it is not produced in the later stages due to its heavy mass. The production of cc 
takes place mainly at the hard vertex. Once produced, cc will evolve in to open D mesons 
and charmonium states, such as J/ip, Xi stc, while travelling through this dense phase of 
matter. Unlike hadronic collisions, where an enhancement of J / is found at tevatron energy 
0, both cc and J/ip get supressed while traveling through different stages of quark-gluon 
plasma. There might be thermal charm productions in the very early stage of the plasma, 
where the temperature is very high, but no production of cc can occur in the later stages 
due to decrease in the temperature. 

The various stages by which the complete evolution of quark-gluon plasma is descrbed in 
URHIC are, i) pre-equilibrium, ii) equilibrium, where one actually studies the thermalised 
quark-gluon plasma, iii) cooling and iv) hadronisation. The production of an equilibrated 
quark-gluon plasma in stage (ii) crucially depends on the pre-equilibrium evolution, i.e. on 
stage (i). Matsui and Satz has suggested the suppression of J/ip in the equilibrium phase. 
In their study, the debye screening length which is calculated from lattice QCD, is found 
to be less than the J/ip radius. This forbids the binding of cc to J/ip in the equilibrated 
quark-gluon plasma. However, such a study is not available for an equilibrating quark-gluon 
plasma. There could be any interesting effects in the pre-equilibrium stage, when cc or 
J/ip travel before reaching the equilibrium stage. As high energy deconfined partons are 
also present in the pre-equilibrium stage, J/ip is supressed due to its interaction with these 
deconfined partons. Since the J/ip dissociation cross section inside a deconfined partons is 
very different from that inside a hadronic gas (see section II), the study of J/ip suppression 
can then provide us with evidence for color deconfinement in the parton plasma and possibly 
with information on the QCD phase transition. 

Another outcome of the pre-equilibrium study is the prediction of the equilibration time, 
the time at which quark-gluon plasma equilibrates. A detail knowledge of equilibration time 
and J/ip formation time is very crucial in determining the J/tp suppressions in URHIC. 
This is because, a large equilibration time may permit a. J/ip formation from cc pairs before 



plasma equilibrates. Then the interaction of a fully formed J/ip with deconfined partons is 
prefered than screening. If the J/ip formation time is greater than the equilibration time 
or comparable to it, then cc interacts with the deconfined partons before screening starts 
operating in the equilibrium stage. The interaction of a cc with the partons before plasma 
equilibriates is not studied yet, and such a calculation is beyond the scope of this paper. 

In any case, a careful study of pre-equilibrium evolution of quark-gluon plasma is the 
consistent way to study production and equilibration of quark-gluon plasma in URHIC. From 
this point of view it is necessary to study what happens actually to the J / ip suppressions at 
different stages of evolution of QGP, rather than estimating it in an equilibrated quark-gluon 
plasma. We study here the suppression of a fully formed J/ if) in equilibrating quark-gluon 
plasma using short distance QCD and suppression of a nascent J/ip by nuclear absorption. 
The J /ip may also be supressed in later stages due to inelastic collisions with hadrons, but we 
do not consider it here. This suppressions will be less due to lesser value of the hadron- J/?/' 
inelastic cross section (described in section II). 

However, J / ip suppression is also measured in p-A collisions where there is no existence 



of quark-gluon plasma |10,11|. At these experiments the suppression of J/ijj is due to 
the presence of a nuclear medium. The prominent mechanism by which these datas are 
explained at SPS is via nuclear absorption, the suppression of a nascent J/ip before it forms 
a physical resonance [T^JT^ . Within this mechanism, one does not assume the existence of a 
quark-gluon plasma phase or a highly densed deconfined partonic medium. This is the case 
with some light nuclei collisions at SPS |]lO|,|lTl. However the recent measurements of J/ip 
suppression by NA50 collaboration |[T^ (Pb-Pb collisions at ^/s = 17GeV) yeild an excess 
suppression of J/ip, which is not explained by the above conventional approach. Presence of 
a nuclear medium may not be enough to explain this data. There are speculations that this 
suppression is due to the existence of a deconfined partonic medium [|15| or a high energy 
density matter |]TB 



Whether an equlibrated plasma has formed in Pb-Pb collisons at SPS(y^ = 17GeV) or 
not, is debatable, but at RHIC(Au-Au collisions at ^/s = 200GeV) we may go closer to the 
equilibrated quark-gluon plasma. In this paper we estimate the J/ip suppressions at RHIC 
and its possible implication on deconfinement. 

However before studying charmonium suppressions in quark-gluon plasma it is essen- 
tial to study the corresponding situations in hadronic collisions. The experimental analysis 
of quarkonia production in high energy hadronic collisions at tevatron energy P| revealed 



a drastic disagreement with the predictions from color singlet model |T^. A substantial 
fraction of quarkonium production at very high energy arises from intermediate color octet 
fluctuation during the evolution into final state color singlet QQ pair. Now a rigorous the- 
oretical foundation has been achived by Bodwin, Braaten, and Leepage [|1^], who devloped 
a QCD formalism by marrying pQCD and an effective field theory within non-relativistic 
QCD(NRQCD). This formalism accounts for the production of both color singlet and color 
octet cc states, that evolve in to final color singlet quarkonium. Within NRQCD, the pro- 
duction amplitude is expanded in powers of both the strong coupling and the velocity 
of heavy quark(f^ = .23 for cc system), which includes higher fock-state components. The 
long distance non-pertubative matrix element is either fitted from experiments or taken, 
in principle, from lattice QCD ||19|. This technique has been used to study many heavy 
quarkonium production processes in hadronic collisons. However a space time evolution of 



quarkonium using this technique has not been achived. This is crucial for studying the 
evolution of J/ ip in nuclear collisions. The existence of background chromoelectric field ||20 



and high dense partonic medium complicates the study of the evolution of cc to J/V', both 
at RHIC and LHC As cc takes ~ 1.0 fm (time in cc rest frame) to form a J/ip bound 
state, it interacts with the nearby light partons and with the background chromoelectric 
field during this period. No attempts is made on this line to study J/ if) production and 
its subsequent evolution in nucleus-nucleus collisions. This might be the consistent way to 
study J/il) suppression in URHIC. 

Since we study the suppression of J / ip in equilibrating quark-gluon plasma, it is relevant 
to discuss the process of equilibration of QGP in URHIC. The equilibration of QGP in 
heavy- ion collisions is not completely studied yet. The rate of equilibration is different for 
diiferent models. Some relevant models that describe the equilibration of QGP are parton 
cascade model(PCM) [^, heavy-ion jet interaction generator(HIJING) |^ and color flux- 
tube model |2^. PCM and HIJING are both perturbative QCD based models, and color 
flux-tube model is a field theoritical model. The former models describe the evolution of 
hard and semihard partons and the later describes the production of soft partons via non- 
perturbative Schwinger mechanism. For a complete study of production and equilibration of 
quark-gluon plasma, color fiux-tube model may be combined with the pQCD based models. 
As an initial attempt, we use here the distribution of partons from parton cascade model to 
study J/ip suppressions in a thermally equilibrating quark-gluon plasma at RHIC. Within 
parton cascade model, a relativistic transport equation is solved to study the distribution of 
partons by considering the direct collision processes 2— > 2, the inelastic collision processes 
2— > 1, and the decay processes 1— > 2 [^. More explicitely, this is written as 

p''dji{x,p) = a{x,p) (1) 

where fi{x,p) is the distribution of a particular type of parton i in the usual phase space. 
The collisions term Ci{x,p) is derived from pQCD taking all the processes mentioned above. 
The initial condition on fi{x,p) is the distribution of partons inside the nucleus before the 



two nuclei collide with each other |2J]. However for a complete study, the addition of soft 
partons and the effect of background field has to be taken into account. These backgruond 
electric field may not exit in the very early stage where pQCD is applicable |2^, but in 



later times this will have an important role in the equilibration of the plasma |2^.2^. Such 



calculations will be available once the color flux-tube model is combined to the pQCD based 



models p3| . The addition of soft partons to hard and semi hard partons within color fiux- 
tube model will enhance the J/ip suppression (see below). This will be taken up separately. 

In section II we describe the short-distance QCD and the dissociation of J/ip. Section 
III contains a brief discussion of nuclear absorpton mechanism. Results and discussions are 
presented in section IV, for RHIC energy. We summarise and conclude the main results in 
section V. 



II. J/ij DISSOCIATION BY DECONFINED GLUONS 



In the framework developed by Bhanot and Peskin p4p8| , interactions between light 



hadrons and deeply bound, heavy quarkonium states, such as the J/ip, are mediated by short- 
range color dipole interactions. For sufficiently heavy quarks, the dissociation of quarkonium 
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states by interaction with light hadrons is fully accounted for by short- distance QCD P7| , p5| . 
Because of its small size, a heavy quarkonium can probe the short distance properties of 
light hadron. A parton based calculation of the J/ip-hadron cross section is thus possible 
via an operator product expansion method, simillar to that used in deeply inelastic lepton- 
hadron scatterings |p7|- p0|] . These pertubative calculations become valid when the space 
and time scale associated with the quarkonium state, tq and tq, are small in comparison 
to the nonperturbative scale Ag^^,, which is the characteristic size of the light hadrons, i.e. 
Tq «< Xqcdj <« '^QCD- For J/tp ground state rj/^ ~ 0.2/m = {lGeV)~^ and 

Ej/^{2Md — Mi/,/) ~ 0.64 GeV. With Xqcd = 0.2GeV, the above inequalities seem to be 
well satisfied and one expect that the dissociation of J/ip by hadron will be governed by 
the J/'?/'-hadron break-up cross section as calculated in short distance QCD. The operator 
product expansion allows one to express the hadron- J/-?/' inelastic cross section in terms of 
the convolution of the gluon- J/-?/' dissociation cross section with the gluon distribution inside 
the hadron. The gluon- J/-?/' dissociation cross section is given by ||27| , |28|j30[| : 



Here Qs is the coupling between gluon and charm quark, rriQ the charm quark mass, and g° 
the gluon energy in the J/ijj rest frame. The heavy quarkonium is coulomb like, and one 
writes, 6o = (|g)2mg PJH. 

The gluons which are soft inside a pion are not capable of dissociating a charmonium. 
Hence they give a low TT-J/ip cross-section. However, deconfined gluons, which carry enough 



energies, are sufficient to break a charmonium [31|. These high energy partons are present 
in URHIC, as observed in several calculations, such as pO| , ^ . These deconfined gluons can 
break a fully formed J/ip that exit inside QGP. This conclusion does not seem to be affected 
substantially by nonperturbative effects ||32l . 



We use the above g-J/i^ dissociation cross section to study the J/i/j suppressions in a 
thermally equilibrating quark-gluon plasma at RHIC. For the central collisions, the J/ip 
survival probability is: 

Sg-j/^iPr) = exp[- J drngi^r) < Vreicr{k ■ v) >k\. (3) 

Here ng{T) is the gluon number density which evolves as system expands longitudinally. 
According to Bjorken scenerio this number density is a function of the boost invariant 
parameter r (r = ^ [t? — z'^)). The thermal average g-J /i^ cross section, < Vreic^ik ■ v) >k, 
is 

Jd'kVreia{k-v)f{k',T{T)) 

< ■ Id^kfikO,T{r)) ■ 

Here v{= {Mt, Pt,0)/Mj/^) is the four- velocity of J/ip in central rapidity region and k 
is the four-momentum of gluon in the parton gas. In the thermal average cross-section 
we use the distribution function, /(A;°,T(r)) = cxp(fcO/T(T))-i gluon. Here a(r) cap- 
tures the deviation from equilibrium, which in principle, is determined from the relation, 
rig(r) = / dT{p^u^)f{p, T{t)), by knowing the evolution of temperature and Ug, for an equi- 
librating quark-gluon plasma. Here dT = (^^Sj'^p "^i^^ 7 = 2 x 8 (the product of spin and color 



degeneracy), and u( = (t/r, 0, 0, z/r)) is the flow velocity. However, as a(r) cancels both 
from numerator and denominator of equation(4), we need not know its form. What mat- 
ters is and T(r) whose evolutions are taken from PCM ||2l]]. Within PCM, %('r) ~ 
0.72n(r), with n(r) = 565/m-3(^)-0-90 and T(r) = 950Me1/(^)-°-30, for a thermally 
equilibrating quark-gluon plasma at RHIC. Here tq = .05/m. Using this distributions the 
integration in equation (3) is performed numerically to study the J/ip survival probability. 



III. J/tp SUPPRESSION BY NUCLEAR ABSORPTION 

Consider p-A collisions: unlike hadronic collisions, where an enhancement of J/ip is 
measured at very high energy , J/-?/; suppression is measured in p-A collisions. A number 
of explanations have been given [|12],|13|,^^ for this suppression, among which nuclear 
absorpton is the prominent. This nuclear absorption is well described according to Glauber 
theory [|12|. Within this theory a J/ip survival probability is written as 



1 f 

SpA = d'^bdzpAih, z) exp(- j dz'pAih, z')aabs{'^ - IM))- (5) 

Here b is the impact parameter where the proton collides with the nucleus and z is the 
longitudinal distance at which a cc pair is created. After the cc is produced, it travels in the 
forward direction colliding inelastically with the surrounding nucleons. The inelastic cc-N 
absorption cross section, aabs, is determined from the experiment. The nuclear density, pa 
at a point r(b, z) is normalised such that: / (i^rp^(r) = A. After integrating over z, we get 

SpA = (-ry^) / d%[l - exp(-o-„b.(l - iM)r^(b))], (6) 

where Ta(s) = pA{SiZ)dz^ is the nucleon density per unit area in the transverse 
plane (transverse to the collision axis), and A' = A{1 — 1/A). Extendng to nucleus-nucleus 



collisions [|T2| we write the survival probability as. 

Sab = J d% j dh [1 - exp(-a,fe,(l - iM)r^(s))] 

X [1 - exp{-aabsil - 1/B)Tb{s - b))] . (7) 

Here b is the impact parameter of collision and s describes the transverse coordinates of 
the interacting nucleons which produce cc. This equation can be written in a more simpler 
form. 

Sab = exp{-aabsP{LA + Lb)), (8) 

where La+Lb is the effective length a cc travels inside two nuclei. Except for the case 
of Pb-Pb collision, where an anomaluos J/ip suppression is observed by NA50 [Il4| , p6| , this 
method explains almost all datas upto S — U collisions, using a common absorption cross 
section ~ 6.2mb. 

For a fixed impact parameter b, this survival probability is; 



Sc-c-Nib) = — [ £s [1 - expi-aaU^ - iM)T^(s))] 

J^A'B'{0)aabs J 

X [1 - exp(-a„,,(l - 1/5)Tb(s - b))] , (9) 

where T4b(6) = J dhTA{s)TB{s - b). 

In our calculation we use a uniform density(i?A = tqA^/^, with tq = 1.2/m), which is 
applicable for heavy nuclei and have adopted a absorption cross section of 6.2mb(used by 
the NA50 collaboration [^). For central Au-Au collisions(6 = 0), we obtain from the above 
equation, Scc-n — 0.42. 



IV. RESULTS AND DISCUSSIONS 

In fig-1 we present the thermal averaged gluon- J/ip cross section as a function of J/V" 
transverse momentum P^, for different values of temperature. It can be seen that this value 
is decreased as the temperature of the plasma is increased. Also for J/ip with high this 
thermal cross section becomes less. This has a direct impact on J/ip suppression. As Pt 
becomes large the survival probability must increase. However, due to large initial parton 
density formed in the initial stage of the plasma, these smaller thermal cross sections do 
not have a large effects on the total survival probability. This is shown in fig-2 where we 
have plotted the survival probability Sg^j/^ as a function of Pt for a thermally equilibrating 
QGP at RHIC within PCM. As can be seen, the suppression is almost 100 percent for J/ip 
with low transverse momentum. Even for J/ip with high Pt{— lOGeV) the suppression is 
about 95 percent. This huge suppression is due to the presence of a highly densed deconfined 
partonic medium at RHIC. As can be seen from equation(3) and fig-1 this suppression will be 
more for a plasma with larger number density and lower temperature. A denser and cooler 



plasma is obtained once the soft partons are added to the hard and semi hard partons p3 
The measurement of such huge suppression will reveal the existence of such a deconfined 
partonic medium, possibly a thermalised quark-gluon plasma, at RHIC. 

However, what one actually measures experimentally is the total survival probability. 
This implies that other sources of J/ip suppression, without the existence of a QGP phase, 
have to be singled out from the total survival probability in order to see the suppressions 
only from this deconfined partonic medium. In secton HI, we have estimated a large sup- 
pression (about 60 percent) due to the nuclear abosorption before a. J /ip \s formed. If the 
estimated suppression due to all sources without QGP phase is much less than the experi- 
mental findings, one may then hope that the rest of the suppressions is due to the existence 
of a highly densed deconfined partons. The experimental measurement of J/ ip at RHIC will 
be able to reveal all the possibilities and may explain the existence of such a deconfined 
partonic medium, possibly a thermalised quark-gluon plasma. 

On the other hand J/ijj enhancement can occur due to an increase in the (thermal) charms 
produced at a very high temperature, especially in the initial stage of the evolution of the 
system p7|j38|] . We recall that screening of J/ip was originally |Q proposed without taking 
production of thermal charms into account, which was argued to be negligible at T ^ 200 - 
300 MeV. However, the system may approach a temperature which is as large as 900 MeV 
initially, and can therefore produce a large number of thermal charms at RHIC. So to make 
resonable estimates at RHIC and LHC, one has to take into account both thermal and hard 
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J/ip. For simplicity, we have not considered here those J/ip which are created from thermal 
charms. 



V. CONCLUSION 

We have found a reasonable suppression (around 60 percent) before a J/ip is fully formed, 
due to the nuclear absorptoin alone and an additional suppression of about 95-100 percent 
due to a deconfined phase. These two mechanisms give a huge suppression, almost 100 
percent, to those J/ip which are produced from primary charm quarks. In principle, the 
total survival probability is written as 

St = Scc-NSg-J/^Sother (10) 

where first two suppressions {Scc-n, Sg-j/^) are due to pre-resonance nuclear absorption 
and g-J/i' dissociation, and the last oiae{Sother) is due to any other possible sources of 
suppression [^^ -^ such as final state comover scatterings. As the hadion-J/ip cross sec- 



tion is smaller than g-J/ip cross section we expect a lesser suppression of J/ip in the 



hadronisation stage. The prominent suppression seems to be due to the nuclear absorption 
and J/tp dissociation in the deconfined partonic system. Hence the measurement of a huge 
J/ip suppression at RHIC may suggest the existence of a partonic plasma. However, many 
refinements, such as space time evolution of quarkonium production in nucleus-nucleus col- 
lision using NRQCD has to be worked out, before unambigous conclusions can be drawn 
for the pre-resonance suppression of J/ip. In any case the contribution from thermal J/ip, 
which may not be neglible at RHIC and LHC, has to be taken into account. It is only after 
incorporating the above features the certainty of screening can be justified. 
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Figure captions 

FIG. 1. The thermal averaged gluon-J/ip cross section < Vreic > as a function of transverse 
momentum Pt at different temperatures. SoUd hne refers to T= 0.2 GeV, upper broken hue 
corresponds to T= 0.4 GeV and lower broken line refers to T=0.8 GeV. 
FIG. 2. The survival probability of J/ip in a. thermally equilibrating plasma at RHIC. 
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